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1. Introduction 

One candidate for Beyond Standard Model physics is the possibility of new strong interactions, 
which might be relevant as an alternative to the standard Higgs mechanism. There have been 
many recent lattice studies of such gauge theories, with larger numbers of flavors than in QCD, 
or alternative fermion representations. For any given theory, one wants to know if the infrared 
behavior is like that of QCD, with a dynamically generated mass gap, or if the theory is conformal 
in the chiral limit. There are many lattice observables which can be used to investigate conformal 
versus QCD-like behavior, among them: the mass spectrum; the flow of the renormalized gauge 
coupling; the presence or absence of finite-temperature phase transitions; the eigenvalues of the 
Dirac operator; and the renormalization group flow of the bare parameters. 

Measurement of the mass spectrum requires large lattice volumes, which becomes particularly 
computationally expensive as one pushes towards to the chiral limit. Further useful information 
can be extracted from these gauge ensembles by measuring the static fermion potential V(r), and 
hence the force F(r) = dV /dr. If a given theory is conformal, in the infrared the potential should 
be purely Coulomb-like. In particular, one can use the QQ scheme, where the renormalized gauge 
coupling is defined via cc qq (r) = r 2 F(r)/CR, where Cr is a representation-dependent constant and 
the RG scale is the fermion separation r. In a conformal theory, the purely Coulomb-like potential 
V(r) in the infrared corresponds to the gauge coupling ct qq (r) flowing to an infrared fixed point. 
Alternatively, if the theory is QCD-like, no such infrared fixed point should appear. 

We have studied at finite lattice spacing the mass spectrum of two particularly interesting BSM 
theories, both with SU (3) gauge symmetry: Nf = 2 fermions in the 2-index symmetric (sextet) 
representation and Nf = 12 fermions in the fundamental representation. For both models, we find 
it difficult to interpret the mass spectrum and related quantities as having conformal behavior with 
a universal critical exponent [1, 2, 3, 4]. By studying in addition the static potential, we test if 
different lattice observables show consistent non-conformal behavior. 

2. Computational details 
2.1 Method 

For all simulations described here we have used the tree-level Symanzik-improved gauge ac- 
tion, with the conventional lattice gauge coupling jS = 6/g 2 as the overall factor in front of the 
Symanzik lattice action. The link variables in the staggered fermion matrix were exponentially 
smeared with two stout steps [5]. The precise definition of the staggered stout action is given 
in [6]. The RHMC and HMC algorithms were used in all runs. For the molecular dynamics time 
evolution we used multiple time scales [7] and the Omelyan integrator [8]. For our mass spectrum 
investigation, we generated ensembles for a large set of lattice volumes and fermion masses, of 
which only a subset is used here. For the Nf = 12 fundamental study, we measure the potential 
V(r) on lattice volumes 48 3 x 96 and 40 3 x 80 at one bare coupling j3 = 2.2 and fermion masses 
ma = 0.01,0.015,0.02,0.025, corresponding to pion masses m % a ranging from roughly 0.16 to 
0.31. For the Nf = 2 sextet work, the lattice volumes are 48 3 x 96 and 32 3 x 64 at /3 = 3.2 and 
ma = 0.003,0.005,0.006, corresponding to m n a from 0.14 to 0.19. Sextet ensembles at j8 = 3.25 
are currently being generated, to study the cutoff dependence. 
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Figure 1: (left) The effective mass V(r,t) and the fitted potential V(r) for the Nf = 12 fundamental model at r — 5 
on a 48 3 x 96 volume at fermion mass ma = 0.01. The quality of the correlated fit is X 2 /^dof = 2.0. (right) Volume 
dependence of the potential V(r) with lattice volumes 48 3 x 96 and 40 3 x 80 at mass ma = 0.01. 

We extract the potential V(r) from measurements of the Wilson loops W(r,t). To improve 
the signal, we use a combination of HYP-smearing of the time-like links, which reduces the self- 
energy of the fermion-antifermion operator, and various levels of 3-dimensional APE-smearing of 
the space-like links, to create a set of operators from which a correlator matrix can be built as input 
to solve a generalized eigenvalue equation [9]. Here we will only show results for one element 
from the diagonal of the conelation matrix. From the Wilson loops, we extract the effective mass 
V(r,t) = —\nW{r,t + l)/W(r,t), which at sufficiently large times can be fitted to a constant V(r), 
including the covariance matrix for the correlated data. We use the double (nested) jackknife 
method for the error analysis: for every outer jackknife sample of Wilson loops, we have an inner 
jackknife loop to determine the covariance matrix for the effective masses, which is included in the 
fit [10]. We bin the data set until we typically have a total of between 10 and 20 bins. In Figure 1 
we show a typical result for the fitting of the effective mass and the determination of V(r). 

2.2 N f = 12 fundamental 

We first discuss the results for the Nf = 12 fundamental model. At the lightest fermion masses 
ma = 0.01 and 0.015, we have simulations on both 48 3 x 96 and 40 3 x 80 lattice volumes. As shown 
in Figure 1, the potential has no visible volume dependence at the smallest mass, hence 40 3 x 80 
is already sufficient to reach the infinite-volume limit at the heavier masses 0.02 and 0.025. This 
gives the potential V(r) in infinite volume at four separate fermion masses. For each mass, we 
parametrize the potential in some form, then study the mass-dependence of the parameters. A 
standard parametrization of the potential is 

V(r)=V -®+Or, (2.1) 

where one includes both Coulomb-like behavior at short distance and string-like behavior at larger 
separation. Because this is a dynamical fermion simulation, at sufficiently large separation the 
string can break and there is no asymptotic string tension. An alternative parametrization is to 
exclude the short-distance data and fit the potential to the form 

V(r)=V + Or (2.2) 
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for larger separations only. We study both parametrizations of the data. 
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Figure 2: Fits of V(r) with and without a/r: (left) on 48 3 x 96 at ma = 0.01 and (right) on 40 3 x 80 at ma = 0.025. 

In Figure 2 in the left panel we show the potential data for 48 3 x 96 at ma = 0.01. The fit 
including the a/r term is over the range 3 < r < 24, and the linear fit without a/r is over the range 
10 < r < 24. The fitted linear parameter is aa 2 = 0.00348(12) and 0.00639(17) for the two fits 
respectively, with the quality of each fit being X 2 /^doi = 32.8/19 and 3.8/13. Note that we do not 
include the correlation in r of the data in the fit. At larger separation the data shows little curvature, 
and the linear fit appears to describe the data perfectly well, without the detection of Casimir energy 
in string formation. To demonstrate that this is not a feature only at the lightest fermion mass, we 
show in Figure 2 in the right panel the potential for 40 3 x 80 at ma = 0.025 . Including the a/r term, 
fitting over the range 3 < r < 20 gives oa 2 = 0.01088(25) with X 2 /^dof = 16.7/15. The linear fit 
without a/r for 10 < r < 20 correspondingly yields aa 2 = 0.01528(55) and % 2 /N Ao f = 2.1/9. 
Again, the potential V(r) at larger separation r is very well described by purely linear behavior. 
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Figure 3: Fits of the string tension for Nf = 12 fundamental theory: (left) from V(r) fits including a/r and (right) 
from linear V(r) fits without a/r. 

Given the determination of a at each fermion mass ma, we now examine the behavior in the 
chiral limit ma — > 0. In a conformal theory which has been deformed by a small fermion mass 
m, quantities with mass dimension, such as particle masses and a 1 / 2 , have a power-like behavior 
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oc m l l y,, \y m = 1 + 7, where / is the anomalous dimension [11, 12]. The critical exponent is univer- 
sal for all particle quantum numbers, and all mass gaps vanish in the chiral limit. Alternatively, if a 
given theory is like QCD with spontaneously broken chiral symmetry, only the Goldstone bosons 
are massless in the chiral limit, all other states are massive, and a 1 / 2 should be non-zero in the 
chiral limit. In Figure 3 we show fits of the mass dependence of a 1 / 2 , testing for QCD-like be- 
havior (parametrized with a linear mass dependence) or conformal power-law behavior. The left 
panel shows fits where a was determined from V(r) fits including the a/r term. Neither linear 
nor power-like behavior describes all four data, hence ma = 0.025 is excluded. Using the three 
smallest masses, the linear fit yields a^ 2 a = 0.0338(23) in the chiral limit, whereas the power-like 
conformal fit gives y = 0.92(12), however both fits are of very poor quality. If instead one uses 
a as determined from linear fits of V(r) at larger r only, the behavior is much improved. Data at 
all four masses can be fitted, and both linear and power-like ansatze fit the data well. The linear 
fit gives a chiral limit value c l l 2 a = 0.0516(23), the power-like fit gives an anomalous dimension 
7 = 1.17(1 1). However, the conformal fits are in very strong tension with the mass spectrum anal- 
ysis. For example, the pion mass dependence indicates a value y = 0.393(3) for the anomalous 
dimension, while the pion decay constant is best described with y = 0.214(16). Given this large 
violation of universality of the critical exponent, we conclude that the indication from the potential 
is that the Nf = 12 fundamental theory has serious problems with the conformal interpretation. 
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Figure 4: Fits of the potential V(r) for the Nf = 2 sextet theory, with and without the a/r term: (left) 48 3 x 96 at 
ma = 0.003 and (right) 32 3 x 64 at ma = 0.006. 

2.3 Nf = 2 sextet 

We next summarize our results for the Nf = 2 sextet model, where the method and analysis 
are very similar to before. We have fewer large volumes and cannot empirically show that volume- 
dependence of V(r) is negligible. Hence we analyze three ensembles: 48 3 x 96 at ma = 0.003, 
and 32 3 x 64 at ma = 0.005 and 0.006. The corresponding pion masses are approximately m n L = 
6.5,5.6 and 6.2, with L the spatial size, giving some indication that volume-dependence should 
be small, as seen in our sextet spectroscopy. In Figure 4 we show fits of V(r), with and without 
the a/r term, for the smallest and largest mass considered. On the largest volume, the data at 
larger separation are again well described by purely linear behavior. In Figure 5 we show fits of the 
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Figure 5: Fits of the string tension for Nf = 2 sextet theory: (left) from V(r) fits including a/r and (right) from linear 
V(r) fits without a/r. In the right plot, the fitted conformal exponent 1/(1 + 7) is consistent with zero, hence the curve 
is omitted. 

mass-dependence of a 1//2 , using both linear and conformal power-like m l ^ l+r ^ forms. As before, 
we consider both parametrizations of V(r). We see in fact very little mass dependence. With 
or without the Coulomb term when extracting a, the fitted conformal exponents are respectively 
1/(1 + 7) = 0.04(4) and 0.00(6), giving unacceptable values of the anomalous dimension 7. (A 
negative value for the exponent would be unphysical and simply reflects statistical fluctuations.) 
Linear extrapolations give a clear non-zero value for the string tension in the chiral limit. This 
suggests that the sextet theory appears to be non-conformal, which is consistent with our analysis 
of the mass spectrum. 

2.4 Force 

In fitting the potential V(r), correlation between data at different r was not taken into account, 
given the instability of the covariance matrix without very large statistics. This can be partially 
cured by extracting the force F(r) directly from the Wilson loops W(r,t). We construct an effective 
force F{f ,t) = V(r + l,t) — V(r,t), which is fitted at sufficiently large time t to a constant. In the 
fit, the covariance matrix includes correlation of the data both in r and in t. The naive definition 
of the force location is r 1 = r + 1/2, which we improve by taking into account the propagator for 
the improved action. For example, in our action r = 4 corresponds to r' = 4.45787, at larger r the 
deviation from half-integer quickly vanishes. If a given theory is conformal, at large r the force 
should have a pure 1/r 2 behavior, such that the renormalized coupling ct qq (r) = r 2 F(r)/CF flows 
to an infrared fixed point with increasing r. Alternatively, linear behavior in the potential V(r) at 
intermediate separation corresponds to a constant force F(r). 

In Figure 6 we show the force as extracted from the largest volume at the lightest mass for both 
the Nf = 2 sextet and Nf = 12 fundamental theories (we find similar behavior at larger mass). As 
the separation r' increases, the force appears to flow to a constant, consistent with the independently 
determined value of a from the potential V(r). We compare with perturbation theory, starting the 
RG flow of a qq from its directly measured value at r' = 3.42522. The perturbative prediction of 
a quickly decreasing force is not supported by the data, and the renormalized coupling continues 
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fundamental N f = 12, P = 2.2 sextet N f = 2, P = 3.20 




Figure 6: The force F(/) compared to 1-loop perturbation theory and a as determined from fitting V(r) (left) Nf = 12 
fundamental on 48 3 x 96 at ma = 0.01 and (right) AT/ = 2 sextet on 48 3 x 96 at ma = 0.003. 

increasing without any indication of an infrared fixed point. The effect of the finite fermion mass, 
and whether or not this behavior is altered in the chiral limit, remains to be explored in future work. 
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